Level Density of a Bose Gas and Extreme Value Statistics 
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We establish a connection between the level density of a gas of non-interacting bosons and the theory of 
extreme value statistics. Depending on the exponent that characterizes the growth of the underlying single- 
particle spectrum, we show that at a given excitation energy the limiting distribution function for the number 
of excited particles follows the three universal distribution laws of extreme value statistics, namely Gumbel, 
Weibull and Frechet. Implications of this result, as well as general properties of the level density at different 
energies, are discussed. 
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The level density is an essential quantity in determining 
the thermodynamic properties of closed quantum systems. In 
interacting many-body (MB) systems its computation is in 
general a difficult problem. The most common framework 
is a mean-field approximation, where a gas of independent 
(quasi-)particles moves in an average self-consistent poten- 
tial. In this case, the energy of the gas is expressed as the 
sum of the occupied single-particle (SP) energies. The com- 
putation of the MB level density thus reduces to a combina- 
torial problem: counting the number of ways into which the 
energy can be distributed among the particles. The level den- 
sity has been extensively studied in fermionic systems, where 
detailed experimental data exists at different excitation ener- 
gies and quantum numbers (see, for some recent progresses 
in this field, Refs.fjJEI])- In spite of the experimental break- 
throughs of the 90's and of the many interesting developments 
that followed, the case of bosonic systems is much less known. 
Studies of the spectral properties have concentrated on the low 
energy range of the spectrum of the condensate phase, where 
collective effects and interactions play a crucial role. 

Our aim here is to compute, within an independent-particle 
approximation, the asymptotic properties of the MB level den- 
sity p MB {E, N) of a Bose gas as a function of the excitation en- 
ergy E and the particle number N. We will first consider two 
extreme regimes that correspond to the degenerate quantum 
and the classical limits of the gas. The level density in these 
two extreme cases behaves quite differently as a function of 
energy. In the former case, where one takes the N — > oo limit 
first keeping the energy E finite, the level density p MB (E, oo) 
increases with energy in a stretched-exponential manner for 
large E. In contrast, in the classical limit where one keeps N 
finite and takes the large E limit, the level density increases 
with energy in a power-law fashion. This leads to a natural 
question: what happens in between these two extreme regimes 
where both E and N are large but finite? The main result of 
this Letter is to show that in this intermediate regime the level 
density displays a rich variety of scaling behaviors depending 
on the SP spectrum and has an interesting connection to the 
extreme value statistics (EVS) of independent random vari- 
ables. 

To explore this intermediate regime, we stay close to the 
degenerate gas limit and compute explicitly the effect of a fi- 



nite number of bosons N on the level density. In this regime, 
in a given configuration of excitation energy E, only a frac- 
tion of the particles contribute to E, the rest remain in the 
ground-state. However, the ground-state occupancy and con- 
sequently the number of excited bosons fluctuate among dif- 
ferent configurations belonging to the same excitation energy 
E. These fluctuations may be small or anomalously large de- 
pending on the SP spectrum. To obtain a quantitative estimate 
of these fluctuations, we compute explicitly the distribution of 
the number of excited particles for a fixed (but large) E. We 
will show that the fraction of configurations at excitation en- 
ergy E with N or less excited bosons, among all possible con- 
figurations belonging to the level E, has a limiting distribution 
(when suitably scaled) for large N and large E. Depending on 
the index v that controls the growth of the SP number of states 
(cf Eq. l|9} below), we show that three limiting distributions 
emerge, namely Gumbel, Weibull and Frechet distributions. 
Interestingly, precisely the same three limiting distributions 
characterize the EVS of independent random variables |3], a 
field that has seen a recent resurgence of interests |4]. 

Our work thus provides a link between these two a priori 
unrelated fields, namely the combinatorial problem associated 
with a non-interacting Bose gas and the EVS. We believe that 
this link is of interest in different branches of physics (such 
as in the computation of black hole entropy |5]), mathemat- 
ics and computer science. For instance, it is well known that 
the computation of the level density of a Bose gas in a one- 
dimensional (ID) harmonic potential (equidistant SP spec- 
trum) is directly related to the theory of partitions of an in- 
teger 1 6| [J ■ The theory of partitions has given rise to deep 
and fundamental results in mathematics, some of them related 
to unique developments by their originality and importance 
|0] . Hence our results also provide a link between the number 
partitioning problem and the EVS, generalizing a theorem of 
Erdos and Lehner [9] (see also fTol flTll > which states that the 
number of summands in a random partition of an integer is 
asymptotically distributed with the Gumbel law. 

We consider non-interacting bosons confined by some 
single-particle potential whose energy levels are ej, j — 
0, 1,2,.. .. We set erj = without loss of generality. Each 
configuration {n,j} of the gas is characterized by an excita- 
tion energy E = X^^li n j £ j an( ^ a tota l number of particles 
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N = YljLo n h w h ere rij = 0, 1, 2, . . . is the occupation num- 
ber of the j-th SP level in that configuration. The level density 
at excitation energy E of a gas of N bosons is given by 



J2 11 ^ 



6[N-J2', 



3=0 



(1) 

The number of excited bosons is simply N cx = N — no = 
J^jLi n j- Since tiq > 0, it follows that N cx < N. Thus, if 
one just keeps track of only the excited bosons, it is an easy 
exercise to show that p M b(E, N) in Eq. Q can alternately be 
interpreted as the number of configurations with energy E and 
with N cx < N. Thus, when N — > oo, p MB (E, oo) simply 
counts the total number of configurations at energy E. 

A convenient way to express Eq. Q is by means of an in- 
verse Laplace transform 

-i />a+ioo />b+zoc 

p MB (E,N) = —^ df3 dae s ^\ (2) 

where 

S(a, (3) = —0Cl(a, 0) + (3E — aN (3) 
is the entropy, 

Af(e) 



dc- 



J3c- a _ i 



the grand potential of the gas, and 

N(e) = ( p{e)de 



(4) 



(5) 



the integrated density of states expressed in terms of the SP 
density of states p(e) = ■ 6(e — €j). In Eq. 0, a and b are 
real parameters such that all the poles of the integrand are to 
the left of the integration path. 

A saddle point approximation with respect to the auxiliary 
parameters a and (3 of the integrals in Eq. yields fl2il 



p UB (E,N) = e s ^/ 2TT^/\V(a,f3)\ 



(6) 



where T>(a, (3) is the determinant of the second derivatives 
of S(a,P). The dependence on N and E in Eq. (|6j arises 
from the saddle point conditions that determine implicitly the 
values of a and (3 in terms of N and E 



jV(a,(3) 



e 0e-a _ I 

jpje) 

e 0e- a _ ! 



N 
E. 



(7) 
(8) 



where Af(a, (3) and £ (a, (3) are the particle number and en- 
ergy functions of the gas, respectively. We will work here in 
the leading order approximation p MB (E,N) w e 5 ( Q,/3 ), and 
thus ignore the prefactor in Eq. ©. 

In Eqs.Q and (|8jl all the non-trivial information is con- 
tained in the SP level density p(e). We use here the continu- 
ous approximation, in which the discreteness of the SP energy 



levels Ej is ignored and p(e) is replaced by a smooth function. 
We assume moreover that the high energy growth of the in- 
tegrated density of states is well approximated, on average, 

by 



Af(e) 



(9) 



Here e is an adimensional energy. To recover dimensional 
quantities in the formulas below, all energies must be mul- 
tiplied by some appropriate factor k. The index v is a real 
positive number that can take arbitrary values depending on 
the confining potential. For instance, if the gas is trapped in a 
one-dimensional potential whose energy levels tj = j s , then 
v = 1/s. In contrast, when the confining potential is a D- 
dimensional harmonic oscillator, then v = D, while when 
it is a D-dimensional box (hard wall cavity potential), then 
v = D/2 (and, for instance, k — (V/6Tr 2 ) 2 / 3 (2m/h 2 ) when 
D = 3, where V is the volume of the cavity and m the mass 
of the particle). 

In the approximation l|9} the weight in the sum ^ of the 
SP ground-state is effectively fixed to zero. Under these con- 
ditions, £(a, (3) (and E) represent the excitation energy of the 
gas, measured with respect to the ground-state energy where 
all particles are in the j — state. For Af(a, (3) (and N) a 
problem appears for v > 1, when condensation may happen. 
In this case, Eq. takes into account only the thermal cloud. 
If needed, we will explicitly incorporate the ground-state oc- 
cupancies in the calculations. 

From Eqs.Q and using (|9jl and consequently p(e) — 
vt v ~ Y , the energy and grand potential are simply related by 
£{a, (3) = —v fl(a, (3). The entropy l|3} may thus be written, 
taking into account the condition (|8j, 



S(a,/3) = (1 + l/v)/3E-aN 



(10) 



For any finite N, a is easily seen from Eq. @ to be negative. 
A standard series expansion of the denominator in Eqs.0 and 
(|8j in terms of z = cxp(a), where < z < 1, allows to 
write, in the continuous approximation, the two saddle point 
conditions as 

Af(a,{3) = l) Li v {z) = N , (11) 

£(<*,0) = VT{ ^ X 1) Li u+1 {z) = E , (12) 

where Li v (z) = YlkLi z k /k v is the polylogarithm function 
and r is Euler function. 

Eq. ( II It shows that Af(a, (3) is an increasing function of 
z. Therefore when N increases at a fixed temperature T — 
j3 , z needs to increase to satisfy the equality. As N — > oo, 
z — > 1. In that limit, the energy is easy to obtain and we get 
£(0,(3) = v J °° rfee7(exp(/3e) - 1) = 6 V /I3 V+1 , where 



6 V = vT(v + l)C(v + 1) 



(13) 



(C(z) = Lii(z) is the Riemann zeta function). From Eq. J12L 
we get the following relation between inverse temperature and 
excitation energy, 



(3 = (3 U = [6 V /E]V^ 



(14) 
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Using this expression for (3 and setting a = in Eq. dlQi . 
we get to leading order in a high energy expansion (i.e., large 
energies compared to the spacing between SP energy levels) 



p MB (£,oo) = exp (1 + l/ V ){e v E v ) 1 '^ 



(15) 



For v = 1 this equation reproduces the well known asymp- 
totic result for an equidistant spectrum ej = j (ID harmonic 

potential), p MB (E, oo) = eV 7 ^/ 6 , obtained by Hardy and 
Ramanujan in the partition problem 1 13] . It was generalized 



to arbitrary 1-D potentials ej ex 



x l v (partitions into non- 
integral powers of integers) in 1 1411 . In the present context, 
Eq. dl5> is valid for any system whose average counting func- 
tion behaves (asymptotically) like Eq. (|5Jl (see also Ref. lllln . 
For instance, it holds for a 3-D harmonic potential (v — 3), or 
a 2-D box of arbitrary shape (or billiard), {v = 1). 

Equation dl5l describes the density in the limit of an infinite 
number of particles for a large but finite excitation energy E. 
In the opposite limit, of a large excitation energy at a fixed 
number of particles, the density behaves quite differently. 
This is the Maxwell-Boltzmann limit, where the gas behaves 
classically. From Eq. dl It . keeping N fixed and increasing the 
temperature (e.g., decreasing /3), it follows that z — > to sat- 
isfy the equality. Then Li„(z) w z for any v, and the station- 
ary phase conditions dl li — J 1 2i become N = T(v + X)zj (3 V 
and_E = vT{v+l)z / (3 U+1 . The relation between temperature 
and excitation energy now is 



E = vNT . 



(16) 



This simple equation generalizes, to an arbitrary confining 
potential, the well-known equipartition of energy valid for 
quadratic Hamiltonians. It provides a precise relation between 
a quantum spectral property (the index v) and the partition 
of energy in the classical limit. From the previous form of 
the stationary phase conditions when z — > we also get 
a = log(/3" N/T(v + 1)). Using this relation for a and 
Eq. d!6t for (3 in Eq. fllOi . the many-body level density now 
takes the form 



(hm(E,N) 



T(v + 1) E v 



N» +1 



N 



(17) 



In contrast to Eq. 1.151 . in the classical limit the level density 
has a power-law dependence on the excitation energy (similar 
results in some particular cases were obtained in 1 15], using 
different methods). When E N 1, using Stirling's 
approximation this equation may be written as p M b{E, N) — 

^"iv' 1 ^ ( ^ ' Under this form, this result coincides for 

v = 1 with the result obtained in Ref.|6] for the asymptotic 
behavior of the partition of integers with a maximum number 
of summands (see also |^). 

So far, we have derived two distinct behaviors of the level 
density with excitation energy: a stretched-exponential behav- 
ior in the quantum-degenerate gas limit, and a power-law be- 
havior in the high temperature classical limit. In the classical 
limit, in any typical configuration of energy E all the parti- 
cles of the gas are excited, while in the quantum-degenerate 



case only a finite fraction of the total number of particles con- 
tribute to the excitation energy (the remaining particles are 
in the ground state). To have a better understanding, in the 
latter case, of the distribution of the number of excited par- 
ticles among all the configurations of energy E, and to gain 
some insight about the transition between the two extreme 
regimes, we now compute, starting from the degenerate-gas 
limit z — > 1, finite N corrections. 

We are interested in particular in computing the relative 
density F(E,N) — p MB {E, N)/p WB (E, oo). This quantity 
gives, among all the possible states of energy between E and 
E + dE, the fraction of those whose number of excited par- 
ticles does not exceed N. Interestingly, we find three distinct 
behaviors for F(E, N), depending on the value of v. In terms 
of a suitable rescaled variable x that depends on N, E and v 
(cf below), the fraction F(E, N) behaves as 



v= 1 : F(E, N) = 
0<i/<l : F(E,N) = 

v>\: F(E,N) = 



exp(- exp(-aj)) , 




x < 



_ x -v/(i-»)-\ x>0 



exp 

exp(— |a;| 7 
1 



x < 
x > . 



(18) 
(19) 



(20) 



In the latter case, the index 7 depends on the precise value of 
v (see Eq. (I23» . These three distributions are known as Gum- 
bel, Frechet and Weibull, respectively. They are the three uni- 
versal limit distributions well known in the theory of extreme 
value statistics of uncorrelated random variables |3]. Below 
we outline the main steps in the derivation of Eas.( ll8l — ( I20l > 
(details will be published elsewhere). 

To prove Eq s . (1 1 8 b— d20l> one needs to compute from 
Eqs.lO and O a(E, N) and /3(E, N), and to replace them 
in the expression (I10i . This is done for large but finite values 
of the particle number N, i.e. in the limit z = c~ v — * 1, 
where 77 = —a is a small positive parameter. Mathematically, 
this requires the computation of the leading order behavior, 
when 77 — * + , of the polylogarithmic function Li v {e~ v ). 
This can be achieved either by a direct computation of the inte- 
grals, or by relying on existing results II 611. Once this is done, 
F(E, N) is obtained by dividing the density by p M s(E, 00) in 
Eq. O- 

Case I: v = 1. We find that the appropriate scaling variable 
for the limiting distribution Eq. dl 81 is 



1 



x = p 1 N + \og/3 1 , 



(21) 



where /3j = (n 2 /6E) 1 / 2 was defined in Eq. JT~4b . It follows 
from Eq. dl 81 that the asymptotic value for the typical number 
of excited bosons for states of energy E is fc 1 log^f 1 . In 
the case of an equidistant spectrum ej = j, this result repro- 
duces the one obtained by Erdos and Lehner in Ref. 1 9] for the 
partition problem. 

Case II: < v < 1. From the procedure described above, 
now we obtain for F(E, N) the Frechet distribution, Eq. dl9l >. 
with the rescaled variable given by 



< v < 1 



N 



(22) 
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where c v = [(1 - u)/v]^- u y u [T(l + v)T(l - 
v)] 1 ' 1 ' /el /[1+v) . Note that in Eq. £3 the exponent 
v/{\ — v) is positive in the corresponding range of v. This 
distribution implies that the typical number of excited bosons 
for states of energy E is c l/ E 1 /( 1+l/ * > /2. However, note that 
the distribution is strongly asymmetric, with a power-law 
decay (toward 1) for N much larger than the typical value. 

Case III: v > 1. This case is slightly more complicated 
than the previous ones, because of the presence of a phase 
transition. In contrast with the previous cases, as N increases 
and z — ► 1 in Eq. (lilt at fixed /3, the function Li v {z) tends to 
a finite value. At constant temperature, there is thus a critical 
number N c = T(l + v)C,(v) / j3 v of bosons that can be hosted 
by the thermal cloud, above which a Bose-Einstein condensa- 
tion starts. We find that the relevant variable in this case is not 
N but the difference N — N c . The behavior of the distribution 
is different according to whether N is smaller or larger than 
N c . When N N c , the exponent 7 and the rescaled variable 
x in Eq. d20t depend on the precise value of v. Three different 
regimes are found, summarized as follows 



\<v<2 : 7 



v = 2 : 7 = 2; x 



v-\ lX \vY{v - l)r(2 - z/)] 1 /" 

1 (JV - N e ) (24) 



/3 U (N-N C ) 



(23) 



v > 2 : 7 = 2; x = 



\og[{ft{N c -N)/v)] 

ft /2 {N-N c ) 
[T{v +1)C(^-1)]V2 



(25) 



where (3 U is given in Eq. (114-1 . Finally, for any u > 1 and 
N > N c (that corresponds to x > 0), a macroscopic fraction 
of the particles is in the ground state. These particles do not 
contribute to the excitation energy, and their precise number 
is unimportant. The behavior of the system is thus identical to 
that of the N — > 00 limit, implying F(E, N) = 1 for N > N c 
(or x > 0). This completes the demonstration of the Weibull 
distribution, Eq. d2"0"l . 



The connection to the number partitioning problem be- 
comes evident if one chooses ej = j and E to be a pos- 
itive integer. The relation E = J^'jLi n jj th en corre- 
sponds to partitioning E into non-zero integers and iV cx = 
Y^jLi n j corresponds to the number of terms or summands 
in a given configuration of partition. The ratio F(E, N) = 
p MB (E, N)/p MB (E, 00) then represents the probability that the 
number of summands in a random partition of integer E is 
less than or equal to N. The corresponding limiting Gum- 
bel law for F(E, N) was first proved by Erdos and Lehner 
by rigorous methods |9[. Our results provide a generaliza- 
tion of this theorem to an arbitrary set of summands charac- 
terized by the growth law Eq. (|9ji- The particular case ej = j s 
with s > corresponds to partitioning an integer E into 
sums of s-th powers of non-zero integers. For example, for 
s = 2, the integer 5 can be partitioned into sums of squares 



as 5 = 2 1 



l z = V 



V + V 



1 + 1 . We have shown 



that while for s = 1 we recover the Gumbel law, the limiting 
distribution of F(E, N) is Frechet for s > 1 (or < v < 1) 
and Weibull for s < 1 (or v > 1). 

In conclusion, we have shown that the density of states of a 
system of independent bosons is described in a suitable scal- 
ing limit by the three limiting laws of extreme value theory. 
This result has a universal character since it depends only on 
a single parameter v that governs the large energy asymptotic 
average behavior of the SP energy spectrum (and is indepen- 
dent, for instance, of the fluctuation properties of the SP spec- 
trum). The derivative of the fraction F(E,A) is related to 
the probability density of the number of excited particles. Our 
general results should allow to recover the moments of the dis- 
tribution of the ground state occupation numbers computed in, 
e.g., Refs. |7]. A probabilistic interpretation of our results in 
the light of 1 17] may shed further light on the connection with 
extreme value theory. 

This work was supported by grants ACI Nanoscience 201, 
ANR NT05-2-42103, ANR-05-Nano-008-02 and the IFRAF 
Institute. 
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